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Geometries in Cellular Network Localization 

Javier Schloemann, Harpreet S. Dhillon, and R. Michael Buehrer 


Abstract —In this letter, we present a new metric for character¬ 
izing the geometric conditions encountered in cellular positioning 
based on the angular spread of the base stations (BSs). The 
metric is shown to be closely related to the geometric-dilution-of- 
precision (GDOP), yet has the benefit of being characterizable 
in terms of the network parameters for BS layouts modeled 
according to a Poisson point process (PPP). As an additional 
benefit, the metric is shown to immediately yield a device’s 
probability of being inside or outside the convex hull of the BSs, 
which localization researchers will widely-recognize as being a 
strong indicator of localization performance. 

Index Terms —Localization, geometric dilution of precision 
(GDOP), hearability, stochastic geometry, point process theory. 


I. Introduction 

ITH the proliferation of smart phones enabling location- 
based services as well as increased pressure due to 
federal regulations such as the FCC E911 mandate [1], cellular 
network localization has recently garnered increasing interest. 
In order to gain insight into how network parameters affect 
cellular localization performance, industry has typically relied 
on simulation-based results. This is because it is a standard 
practice to model BS deployments using a regular grid model, 
which does not lend itself well to mathematical analysis. 
Motivated by the current capacity-driven irregular deployments 
of low-power base stations (BSs), such as pico and femtocells, 
there have been several successful attempts at characterizing 
the downlink and uplink performance of cellular networks by 
modeling the BS locations according to a homogeneous Poisson 
point process (PPP) (e.g., [2], [3]). In addition to being more 
realistic in capturing the irregularities of current deployments, 
the PPP-based models lend tractability by allowing the use of 
tools from stochastic geometry in network analyses. 

The PPP approach is the basis of a new model proposed 
in [4], [5] for studying cellular positioning. The three funda¬ 
mental components determining localization performance are 
(i) the number of participating BSs, (ii) the locations of these 
BSs relative to the device being localized, and (iii) the accuracy 
of the positioning observations. In [5], the model was used to 
obtain accurate analytic expressions of BS hearability , i.e. the 
number of participating BSs. A well-known and fundamental 
metric for evaluating the geometric component is the geometric- 
dilution-of-precision (GDOP) [ 6 ]. This metric is closely related 
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to localization performance, but is limited to studying fixed 
positioning scenarios. Instead, it is desirable to understand how 
network parameters affect the favorability of all geometries 
which may be encountered. This motivates the need for a new 
tractable metric which can be used to gain insights into how 
network parameters affect localization geometries. 

II. System Model 

A. Spatial base station layout 

The locations of the BSs are modeled using a homogeneous 
PPP $ G R 2 with density A [7]. Due to the stationarity of a 
homogeneous PPP, the device to be localized is assumed to be 
located at the origin o. For the link from some BS cc £ $ to 
the origin, the received signal-to-interference-plus-noise ratio 
(SINR) is expressed as: 
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where P is the transmit power, S z denotes the independent 
shadowing affecting the signal from BS z to the origin, a > 2 
is the pathloss exponent, and a 2 is the noise variance. The 
BSs are assumed active with probability /, which represents 
the average network load. To model this loading, a y are 
independent indicator random variables modeling the thinning 
of the interference field by taking on value 1 with probability 
/. Furthermore, localization systems often include some sort 
of multiplicative processing gain 7 due to integration, which 
is assumed to average out the effect of small scale fading. 

B. Number of participating base stations 

For localization, it is well-known that including an increasing 
number of BSs in the localization procedure results in a general 
improvement in positioning accuracy. Thus, we assume that 
for purposes of localization, a device will take advantage 
of as many BSs as it can successfully detect (or hear), 
i.e., all BSs whose signals arrive with some minimum link 
quality. Specifically, a BS 1 participates in the localization 
procedure when 7 - SINR T > /?, where /3 is the post-processing 
SINR threshold above which the signals from the BSs must 
arrive in order for them to be successfully detected. Let 
IN = 1 (SINR. t > (3/ 7 ) represent the number of BSs 

hearable at the origin. In this work, we employ the hearability 
results derived in [4], [5] so that 

P(IN = L) = P L (/, f, ol, (3, 7 , A) - P L+ 1 (/,/,a,/3,7, A), (2) 

where A = AE[iSf /, “] is the shadowing-transformed density 
of the original PPP <I> and P L (/. /, a, 3, 7 , A), the probability 
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of hearing at least L BSs in interference-limited networks 
(including cellular networks), is obtained from the hearability 
expressions presented in [5, Theorem 2 and its corollaries]. 


III. Geometric Characterization 


A particularly nice property of the GDOP is that it is closely 
related to the Cramer-Rao lower bound (CRLB), the lower 
bound on the mean square error (MSE) of an unbiased position 
estimator. For TOA-based localization, if the errors in the 
distance measurements between the mobile and different BSs 
exhibit a common standard deviation of ay meters, then the 
CRLB may be calculated from the GDOP as follows [6]: 

CRLB = of • GDOP 2 (3) 


More specifically, for two-dimensional positioning with position 
[x,y] T , we can write GDOP as 
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where of and of are the variances of the estimates of x and y, 
Vi is the distance to BS i. While it is not difficult to calculate 
the GDOP for specific scenarios, we are not aware of any 
technique for characterizing the distribution of GDOP over all 
positioning scenarios which may be encountered for a given 
set of network parameters. 

In order to characterize the geometric conditions of BSs, we 
now introduce an alternate metric to the GDOP which is based 
on the angular spread of the BSs from the mobile device. 


Definition 1 (Maximum angular separation metric). Consider L 
participating BSs, where the angle from the origin to the /' th BS 
is 0£, measured from some common baseline. For convenience, 
we number the BSs in order such that 6>i < d 2 < ■. ■ < Ol and 

let /.£ = 6e+i — 9i for 1 < i < L and Zy = 27r — (8l — 9\). 
The metric 


^max = max{Zi,..., 4 l} (6 ) 

is then the maximum angular separation between the BSs. 


Now, noting that — 
can re-write GDOP as 

GDOP = 
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Thus, for TOA we find that there is clearly a relationship 
between Z max and GDOP. In this letter, we consider downlink 
TDOA positioning, e.g., OTDOA in LTE [8], which is typically 
the first-choice fallback technique in cellular networks when the 


TABLE I 

Correlation Between Z max and GDOP 


Correlation Metric 

L = 4 

L = 5 

L = 6 

L > 4 

Spearman’s p rank 

0.908 

0.927 

0.924 

0.912 

Pearson’s, GDOP 

0.029 

0.513 

0.659 

0.023 

Pearson's, log(GDOP) 

0.843 

0.870 

0.880 

0.871 


Global Positioning System is unavailable. Showing a clean rela¬ 
tionship similar to (7) for TDOA is slightly more complicated 
than TOA. Thus, in order to compare the distributions of the 
two metrics for TDOA, we take an alternate correlation-based 
approach in which we simulated 2 million TDOA positioning 
scenarios with / = 1, a = 4, / 3/7 = — lOdB, log-normal 
shadowing with a standard deviation of oy = 8 dB, and a BS 
density of A = 2/(yfZ ■ 500 2 ) m 2 which is equivalent to that 
of an infinite hexagonal grid with 500m intersite distances 
(ISD ). 1 In Figure 1, the distributions of the GDOP are plotted 
for all scenarios whose Z max values fell within the specified 
ranges. Small values of Z max (approximately 37 r /4 and below) 
are tightly coupled with low GDOP values. As Z max increases, 
the actual GDOP values are less predictable, though it is clear 
that higher Z max values are tied to higher GDOPs (in fact, the 
GDOPs can be lower bounded). An alternate way to look at 
the connection between the metrics is through a correlation 
analysis. In Table I, the Spearman rank correlation [10] and 
the well-known Pearson’s product-moment (linear) correlation 
between Z max and the GDOP are shown for several values 
of L. First, we note that the linear correlation between Z max 
and GDOP is quite poor. This is not unexpected, considering 
that it is difficult to draw a linear relationship between one 
metric whose values are limited to between 0 and 27 t and a 
second metric which is unbounded. On the other hand, we 
note that the Spearman rank correlation is quite good (above 
0.9 for all values of L). This implies several things: (i) that 
Zmax is a good metric for comparing the geometries of two 
positioning scenarios and (ii) that the relationship between 
Zmax and GDOP can be modeled well by some monotonic 
function. Though determining the best monotonic function 
which relates Z max and GDOP is outside the scope of this 
paper, consider the Pearson correlation between Z max and 
log (GDOP) shown in Table I. The fact that there is a significant 
coupling implies that there exist some values of m and b for 
which GDOP = e mZ,,li ’ x+/ ' provides a simple and reasonable 
approximate transformation between the two metrics. 

Now that a relationship between the two metrics has been 
established, we continue our analysis by characterizing the 
distribution of Z max . What we were seeking is the probability 
that the maximum angle made by the base station is greater than 
some value. Calculating this probability from first principles 
is not straightforward. Flowever, (and quite fortunately) this 
probability is strongly related to two classic problems in 
probabilistic geometry, the probability that L random arcs 
cover the circumference of a circle and the random division 

'Although the analysis is valid for all values of a > 2, we present 
results using a = 4 due to its proximity to the a = 3.76 value used in 
3GPP positioning studies [9] and the availability of simplified expressions for 
hearability in [5] for this special case. 
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Fig. 1. Relationship between the GDOP and Z max in cellular TDOA 
positioning. For small values of Z max , the GDOP is very predictable, while 
larger values of Z max provide a good indication of the minimum geometric 
dilution that can be expected. 



Fig. 2. THE IMPACT OF THE NUMBER OF BSs: As more BSs participate in the 
localization procedure, the distributions of Z max become increasingly more 
desirable (i.e., there is a general reduction in the Z max values encountered). 


of a unit interval (see [11] for a discussion). Concerning the 
former problem, Stevens derived the probability that n random 
arcs of size a covered a unit circumference circle in what many 
consider to be an ingenious proof [12]. It was later shown that 
this result could be related to the distribution of the number 
of arcs needed to cover a unit circle and the random division 
of a unit interval (see [13]). Although space doesn’t allow for 
a full derivation, Stevens’ original result can be modified to 
provide the CDF of Z max . Specifically, for a fixed value of L, 


P(Zmax<^|W= J L)=^ (f)(-l) n ( 

n —0 ' ' 


i-m 

2irJ 


L -1 
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where \ = min{L, \_2'k/^\] [11]. In Figure 2, the theoretical 
distribution of P(Z max < <p|IM = L) is plotted for several dif¬ 
ferent values of L and verified to match the Z max distributions 
observed in the aforementioned simulation. As expected, higher 
values of L are desirable for reducing both the minimum and 
maximum values of Z max . Similarly, higher values of L are 
known to yield generally lower GDOP values. 

Since the distribution of the number of hearable BSs depends 
on the network parameters, it is worthwhile to characterize the 
distribution of Z max across different L values weighted by the 
probability of hearing L BSs. Consider only those scenarios 
where DM > 4, i.e., those times when sufficient BSs are hearable 
to ensure unique localizability using TDOA [14], Then, 

P(Z max < ip, N > 4) 


P(Z max < <£>|IN > 4) = 
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P(Z max <^|N =L)P(IN =L) 
P(DM > 4) 


(9) 


Note that although the summation limits are infinite, it is not 
difficult to calculate (9) with a high accuracy since P(BM = 
L) — > 0 for even relatively low values of L (e.g., see Figure 11 
in [5]). While this behavior eases computation, it makes it 


increasingly more difficult to obtain lower values of Z max as 
the expected number of BSs required to satisfy Z max < ip is 

[tp/2ir- 1] / _ rup\n-1 

E[L]= V (-!)"■ 1 y + i , 
n= 1 V 2 ^r ) 

which grows on the order O (< p _1 logtp -1 ) as <p —> 0 [13]. 2 

IV. Numerical Results and Discussion 

In this section, we present some insights into how the 
network parameters affect the distribution of Z max and also 
reveal a particularly valuable geometric insight which is easily 
obtained from the Z max metric. First, let us consider how 
altering the required pre-processing SINR, P/'y, which is a 
measure of detection sensitivity, affects the distribution of Z max . 
In Figure 3, (9) is plotted for target pre-processing SINRs 
between —6dB and — 14dB (the range of values suggested at 
different points throughout 3GPP positioning studies). While 
hearability clearly improves as shown in [5], the effect on the 
distribution of Z max is surprisingly small. In Figure 4, we again 
consider (9), but now versus changes in the average network 
load /. We see that reducing the network load has a significantly 
greater effect on the geometric conditions measured using Z max . 
Consider the values at P(Z max < 7r|DM > 4) highlighted by the 
dashed lines in Figures 3 and 4. As the network load decreases 
from / = 1 to / = 0.1, P(Z max < 7r| DM > 4) increases 
from 0.5 to 0.95. This is significant because by evaluating 
P(Z max < 7r|DM > 4), one immediately obtains the probability 
of being inside the convex hull of the BSs, which is known to 
be a strong indicator of localization performance [15], [16]. 
This is demonstrated in Figure 5 for both L = 4 and L > 4, 
which reveals a drastic increase in the GDOPs when the device 
to be localized is outside the convex hull of the BSs compared 
to when it is inside the convex hull, where the GDOP values 
are essentially upper-bounded by 4. 

2 The term O denotes the standard big-0 notation which is used to describe 
the asymptotic behavior of a function. 
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Fig. 3. The impact of receiver SENSITIVITY: While being able to 
detect BS signals at lower SINRs generally improves hearability, the effect on 
geometry is not as pronounced as might be expected (/ = 0.5, a = 4, a s = 
8dB, A = 2/(V3 ■ 500 2 ) m" 2 ). 



Fig. 5. Localization performance inside and outside the convex 
HULL: This figure illustrates the value of being inside the convex hull, where 
GDOP values are almost certainly below 4 (f3/y = — 10dB,a = 4, cr s = 
8dB, A = 2/(^3 ■ 500 2 ) m" 2 ). 



Z max < radianS > 


Fig. 4. The impact of average network load: Reducing the network 
load yields a significant improvement of the distribution of Z max . In particular, 
by considering Z ma x < 7r, we observe a drastic improvement in the probability 
of being inside the convex hull of the BSs (/3/y = — 10dB,a = 4, cr s = 
8dB, A = 2/(V3 ■ 500 2 ) m" 2 ). 


V. Conclusion 

In this letter, we presented an alternate metric to the well- 
known GDOP for characterizing the geometric conditions of 
the base stations in cellular localization. The metric is based 
on the maximum angular separation between BSs, which was 
shown to be closely related with the GDOP. In contrast to the 
GDOP, which is limited to studying the geometric conditions 
of fixed positioning scenarios, our metric can be tractably 
characterized in terms of the network parameters across all 
positioning scenarios. Additionally, the metric can be employed 
to immediately obtain the probability that a device will be 
inside or outside the convex hull of the BSs, which is a widely- 
recognized indicator of localization performance. 
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